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Motivations 


o New synth sounds > 
o Low computational cost 
o Virtual Analog Oscillators 
o Usage as audio effect 
o The challenge: > 

“When I first got some - I won't call it music - sounds out of a 
computer in 1957, they were pretty horrible. (...) Almost all the 
sequence of samples - the sounds that you produce with a digital 
process - are either uninteresting, or disagreeable, or downright 
painful and dangerous. It’s very hard to find beautiful timbres.” 

Max Mathews, 2010. 
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Contribution 


LTV theory approach to distortion techniques 
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Phase Distortion 

Phaseshaping - US patent 4658691 Casio - CZ 


Add a phase distortion function to the regular phase generator 
Sawtooth: Inflection point on the regular (dashed) index 


t + g{t) 


0.5)), 0 <t<d 

0.5|5^ + 0.5, d<t< 1 


For d = 0.05 
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The allpass filter 


H(z) 



Flat magnitude response 


Frequency dependent phase shift 


. . „ i / — asin (ui) \ 

4>{u) = -u + 2tan -- r ^- 

— a cos (ijj) J 


Reverb, chorus, flanger, phaser, spectral delay 
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Amplitude modulation 


cos ( 2irf c n ) cos (2irf m n) = ^ cos (2nf c n + 2irf m n )+^ cos ( 2irf c n — 2irf m n) 
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Allpass filters coefficient modulation 





Jussi Pekonen, 2008 


Coefficient-modulated first-order allpass filter as distortion effect 
o Suggests the method for sound synthesis and audio effects 
o Recall that classic PD is restricted to cyclic tables 
o Derives stability condition 

|m(n)| <1 Vn 


o Recommends appropriate values for m(n) 
o Allpass Dispersion on low frequencies 


<t>Dc{n) 


1 — m(n ) 
1 + m(n) 
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Spectrally rich phase distortion sound synthesis using allpass filter 


Time-varying allpass transfer function 


H{z , n) 


— m(r) + z 1 
1 — m(n)z” 1 


Phase distortion 


,. . „ i / - m(n ) sin ftu) 

<Hw, n) = —uj + 2 tan -- — - 

— m{n) cos (a;) 
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4>(lj, n) = — u + 2 tan 
Knowing 4>(u, n ), use tan(x) & x, 


— m(n) sin (tu) 

1 — m(n) cos (w) 


m{n ) 


2 sin (a;) — (0(cj, n) + uo) cos (a;) 
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Emulate the classic phase distortion technique 


t + g(t) 


0.5^, 0 <t<d 

0.5|5^ + 0.5, d < t <1 


Subtract linear phase from the phase distortion function 


g(t) 


ik~ d )b °< t < d 

(| - d )jE% + 0.5, d < t < 1 
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Allpass filters coefficient modulation 

J.Timoney, V.Lazzarini, J.Pekonen, V.Valimaki 

Range for the allpass modulation should be [—cu, —7r] 
u \ g(t)((l — 2d)n - lo) 

{) =-(I^2d)t- (1 ” 2d)1 “ 

Get the modulation function 

m tn) =_ -(#*»)+ _ 

2 sin (w) — (<p(uj, n ) + oj) cos (w) 

Implementation with difference equations 

y(n) = x{n - 1) - m(n)(x(n) - y{n - 1)) 

> 
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Allpass filters coefficient modulation 

J.Timoney, V.Lazzarini, J.Pekonen, V.Valimaki 


Phase distortion and coefficient modulation functions 




Conclusions 

O 


13/25 







Conclusions 


Context Review Dynamic PD FBAM 

OOOOOO0OOOO oooooo 

Allpass filters coefficient modulation 

J.Timoney, V.Lazzarini, J.Pekonen, V.Valimaki 


Outputs with classic PD (solid) and modulated allpass (dashed) 
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J.Timoney, V.Lazzarini, J.Pekonen, V.Valimaki 
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Classic PD and Modulated allpass spectra 
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Allpass filters coefficient modulation 





Arbitrary distortion function 


y(n) = 0.4 cos (f Q ) + 0.4 cos (2 f 0 - -J + 


7r 


7T 


0.35cos [3fo + — + 0.3cos 4 fo + — 


47t\ 


3 J 


Shift it to the appropriate range 


y s («) 


Create your own (: 

o 


7T (y(n) + 1) 
2 2 
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Allpass filters coefficient modulation 





Arbitrary distortion function 


Phase distortion and derived modulation functions 



Samples Index 


Allpass Coefficient Modulation Function 
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Allpass filters coefficient modulation 





Arbitrary distortion function 


Waveform and spectrum 



Time (Samples) 


Spectrum of Modulated Allpass Filter Output 



Frequency (Hertz) 
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J.KIeimola, V.Lazzarini, J.Timoney, V.Valimaki, 2009 


FeedBack Amplitude Modulation (FBAM) 

Revisiting of an old idea by A.Layzer 
tested by Risset in the catalogue > 

Modulate oscillator amplitude using its output 
y(n) = cos (w 0 n)[l + (3y(n - 1)] 
with uq = 2nfo and y[0] = 0 
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FeedBack Amplitude Modulation 


y(n) = cos(w 0 /i)[l + y(n - 1)] 
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= Er=onLoC°sN(n-™)] 


cos 2 (p) = |(1 + cos(2p)) 
cos 3 (p) = ^(3 cos(p) + cos(3p )) 
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LPTV interpretation 

cos(coon) 


y(n) = x(n) + /3a(n)y(n — 1) 

x(n) = a(n) = cos (won) 
in this case (but could be /) 



■f—> Out 



1 pole coefficient modulated HR —» Dynamic PD 
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Feedback Amplitude Modulation 


f3 similar to FM’s modulation index 
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Feedback Amplitude Modulation 

Stability condition 

N 

(3 cos (u^o mi) < 1 

m= 1 

Aliasing before instability 
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2nd order FBAM 

Two previous outputs with individual /3s 

y(n) = cos (w 0 n)[ 1 + /3i y(n - 1) + /3 2 y(n - 2)] 
Narrower pulse and wider band > 



Frequency (kHz) 
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Conclusions 


o Reissue of a classic technique 

o Different kind of implementation 

o Input and modulation can be arbitrary signals 


o Deeper investigation of LTV 
o Studying 2nd and higher order systems stability 


Thanks a lot! 
antonio.goulart@usp.br 
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